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GABOR ANALYSIS FOR SCHRODINGER EQUATIONS AND 
PROPAGATION OF SINGULARITIES 


ELENA CORDERO, FABIO NICOLA, AND LUIGI RODINO 

Abstract. We consider the Schrodinger equation 

du 

i——\-Hu = 0, H = a(x,D), 
ot 

where the Hamiltonian a(z), z = (a:, 77 ), is assumed real-valued and smooth, 
with bounded derivatives |i9“a(z)| < for every |Qf| > 2, z G For such 

equation results are known concerning well-posedness of the Cauchy problem for 
initial data in and local representation of the propagator by means 

of Fourier integral operators. 

In the present paper we give a global expression for in terms of Ga¬ 
bor analysis and we deduce boundedness in modulation spaces. Moreover, by 
using time-frequency techniques, we obtain a result of propagation of micro¬ 
singularities for 


1. Introduction 

Time-frequency Analysis, also named Gabor Analysis cf. I221E6], has found im¬ 
portant applications in Signal Processing and related problems in Numerical Anal¬ 
ysis, see for example |81I12] and references therein. More recently, time-frequency 
methods have been applied to the study of the partial differential equations, in par¬ 
ticular constant coefficient wave, Klein-Gordon, parabolic and Schrodinger equa¬ 
tions |2l|3lE2l|2llEIll32lE3l|36lS5lll6], let us also refer to the survey [H] and the 
monograph m- The analysis of variable coefficient Schrodinger type equations 
was carried out in [U |9l UHl CHI HH UHl [161 UHl US] , see also [171 [TH] in the analytic 
category. 

In the present paper we address to the Schrodinger equation 

{ du 

i— + a{x,D)u = 0 
u{0,x) = uo{x). 
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where the symbol a{z), z = is real-valued and smooth, satisfying 

(2) \d‘^a{z)\ < Co,, |a|>2,2eM2d_ 

In (II]) we understand 

a(x,D)f{x)= f e“<"'Oa(i,0/K)rfL 

jRd- 

with 

/( 0 = [ 

JM.d 

As basic example one can consider the case when a{z) is a quadratic form in z, 
including the free particle operator a{x, D) = —A and the harmonic oscillator 
a{x, D) = —A -I- \x\‘^. 

Under the assumption ([2]), it is easy to show by energy methods that the Cauchy 
problem ([1]) is well-posed in and see for example |13] and we may 

consider the propagator 

(3) SiW^), H = a{x, D), 

mapping the initial datum uq to the solution u{t,x) at time t G M. 

Going further, one would like to obtain an explicit expression for from which 
one may obtain precise estimates for the solutions, and numerical analysis. To this 
end, by using Fourier methods, one expects for small values of f a representation 
as type I Fourier integral operator 

( 4 ) {C^^uo){t,x)= [ 

jR<i 

In fact, under our assumption ([2]), the phase function $ has quadratic growth with 
respect to x,r] and the amplitude b is in the Hormander class Sqq, i.e., bounded 
together with its derivatives. The case of Hamiltonians a(x, r]) of polyhomogeneous 
type was hrst treated in [3 ESI EH]. For the present general case, where homogeneity 
is not assumed, see e.g. [DElEiniEaEgiMlESlESlISl]. 

Time-frequency analysis enters the picture at this moment, giving a global ex¬ 
pression for Let us recall the notation for the time-frequency shifts: 

(5) x,r],t e z = {x, rj). 

The short-time Fourier transform(STFT) of a function or distribution / on with 
respect to a Schwartz window function g G \ {0} is defined by 


( 6 ) Vgf{x,r]) = {f,7r{z)g) = f{v)g{v-x)e^^dvx)(^y^ z = {x,r]) e R^"^. 
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The time-frequency representation of a linear continuous operator P : —)■ 

is provided by the (continuous) Gabor matrix 

(7) k{w, z) := {P7r{w)g, 7r{z)g), w,z e 
so that 

(8) Vg{Pf){z)=[ k{w,z)Vgf{w)dw. 

Following the pattern of |20], we study the Gabor matrix k{t, w, z) of the propagator 
Its structure is linked with the Hamiltonian held of a{x,^). Precisely, consider 

{ 27rx = —V^a(x, .^) 

27ri = VM^,0 
3^(0) = y, ^(0) = T], 

(the factor 27 r depends on the normalization of the STFT). The solution Xt{y^ v) = 
{x(t,y,ri),^(t,y,ri)) exists for all f G M thanks to our hypothesis, and dehnes a 
symplectic diffeomorphism xt '■ The components of Xt are functions 

with bounded smooth derivatives of any order in R x R'^. 

Theorem 1.1. Let k{t,w,z) be the Gabor matrix of the Schrodinger propagator 
Then for every s > 0 there exists C = C{t, s) > 0 such that 

(10) \k{t,w,z)\<C{z-Xt{w))~% z = (x,0, = (i/,il) e R^'^. 

For t small enough our assumptions yield det ^{t, y,g) ^ 0 in the expression of 
Xt, and (ITU]) is then equivalent to (|1|) with the phase <F linked to Xt as standard 
and b{t, ■) G S'q q, see the next Section 2. In the classical approach, cf. |l| . the 
occurrence of caustics makes the validity of (|1|) local in time and for global time 
f G R multiple compositions of local representations are used, with unbounded 
number of variables possibly appearing in the expression. Observe instead that 
k{t,w,z) keeps meaning for every f G R, and the estimates flT0|) hold for Xt with 
t G R. 

A natural functional frame to express boundedness and propagation results for 
g®tR jg gjyg]^ py modulation spaces, see [23] and the short survey in Section |2i 
For l<p<oo,rGR, the modulation space Mf(W^) is defined by the space of 
all / G iS'(R'^) for which 

(11) ll/HWo= / \VJi^W{^rdz<oo 

(with obvious changes for p = cxo). 

From Theorem 11.11 it is easy to deduce the following 
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Theorem 1.2. For every r G M, 1 < p < oo, t G M, we have 


( 12 ) 


e' 


,itH . 




The proofs of Theorems 11.11 and 11.21 are given in [20] when the symbol a{x, rj) is 
a polyhomogeneous symbol. The proof for the present non-homogeneous situation 
follows closely the one in [20]. Our aim being to introduce non-expert readers to 
the methods of the time-frequency analysis, we shall reproduce in the sequel the 
main lines of the argument. 

Novelty with respect to 121 will be a result of propagation of (micro) singularities 
for the solutions of ([I]). In fact, in [20] by taking advantage of the homogeneous 
structure, the propagation was expressed in terms of the global (Gabor) wave front 
set, whereas in ([1]), ([2]) homogeneity is lost in general. This will require the use of 
a more rehned notion, namely the filter of the Gabor singularities, see below. 

Let us first recall that the global wave front set was introduced by Hormander 
[30] in 1991, see also [iQ], where the name of Gabor wave front set was given. 
The work of Hormander [30] was addressed to the study of the hyperbolic equa¬ 
tions with double characteristics, and also provided propagation of singularities 
for Schrodinger equations ([I]) with quadratic Hamiltonian a{z). Such result was 
generalized to different classes of linear and nonlinear equations, beside [20] see for 
example [371139] and [18], concerning the analytic category. One may hnd in these 
papers references to the wide previous literature on the subject. 

The renewed and increasing interest for the Gabor wave front set derives from 
the fact that, under the action of a metaplectic operator, it moves according to 
the associated linear symplectic transformation. More generally, if P is a Fourier 
integral operator as in (jl]) with a phase function <F homogeneous of degree 2 in 
{x, rj), then the Gabor wave front set is determined by the corresponding map Xt- 

Returning to the present non-homogeneous context, our definition of Gabor (mi¬ 
cro) singularity will be as follows. Let T be a subset of R^*^. Given a distribution 
/ G iS'(R'^), we say that / is regular in T, 1 < p < CX3, r G R, if there exists a 
neighborhood T^ of T (see precise definitions in Section 4) such that 



We shall call filter of the Mf singularities of f the collection of subsets of R^'^ 
(14) PP(/) = {A C R^^ / is MP regular inT = R^^ \ A}. 


Note that / is regular in any bounded set T C R^'^. 

We may now state our main result. 

Theorem 1.3. For every r G R, 1 < p < oo, t G R, mq ^ iS'(R‘^), we have 


(16) 


xA.r^(u„)) = JF?(e“"ti„). 
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The content of the next sections is the following. Section 2 is devoted to some 
preliminaries, concerning properties of the short-time Fonrier transform, modnla- 
tion spaces and canonical transformations. Section 3 presents the theory of the 
Fonrier integral operators in terms of Gabor analysis, cf [10]. The proofs of The¬ 
orems 11.11 and 11.21 are given as application. Section 4 contains the analysis of the 
Gabor singnlarities and the proof of Theorem 11.31 

2. Preliminaries 

In what follows there are the basic concepts of time-freqnency analysis. For 
details we refer to [26]. We also discnss the properties of phase fnnctions and 
canonical transformations. 

2.1. The Short-time Fourier Transform. Given a distribntion / G and 

a Schwartz fnnction g G \ {0} (the so-called window), the short-time Fonrier 

transform (STFT) of / with respect to g is defined as in ([6]). The short-time 
Fonrier transform is well-dehned whenever the bracket (•, •) makes sense for dnal 
pairs of fnnction or distribntion spaces, in particnlar for / G and g G 

f,gE Gonsider f,gE then Vgf G 

We recall the following pointwise ineqnality for the short-time Fonrier transform 
[2S1 Lemma 11.3.3], used to change window functions. 

Lemma 2.1. If go,gi,'y G such that {'y,gi) ^ 0 and f G then the 

inequality 

\v,J(x,()\ < |y;W(IL./l»IL.7l)(i,?) 

holds pointwise for all {x,f) G 

2.2. Modulation spaces. Weighted modulation spaces measure the decay of the 
STFT on the time-frequency (phase space) plane. They were introduced by Fe- 
ichtinger in the 80’s [23]. 

Weight Functions. A weight function v is submultiplicative if v{zi + Z 2 ) < 
v{zi)v{z 2 ), for all Zi,Z 2 G We shall work with the weight functions 

(16) y^(2;) = (^)* = (1 + |2;|2)f^ S G R, 

which are submultiplicative for s > 0. 

For s > 0, we denote by A4u^(R^^) the space of Us-moderate weights on R^*^. 
These are measurable positive functions m satisfying 

m{z + w) < Cva{z)m{w) 

for every z,w & R^'^. 
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Definition 2.2. Given g G s > 0, m G and I < p,q < oo, the 

modulation space consists of all tempered distributions f G such 

that Vgf G (weighted mixed-norm spaces). The norm on MfffEf) is 

(17) ll/IU- = rjiu- =( [ ( [ IVgfix^Ol^M^.O^dxY^'dA 

yjR’i \jR>i j ) 

(obvious modifications if p = oo or q = oo). 

When p = q, we write instead of and for M^^(R'^). 

The spaces are Banach spaces and every nonzero g G Mj|^(R'^) yields 

an equivalent norm in (ITTD and hence is independent on the choice of 

g G MlfiMfi). 

We recover the Hormander class 

(18) (*''')■ 

s>0 

For any 1 < p, g < cx), 


(19) Pi = 5(R‘^), IJ MP:‘((R‘^) = 5'(R'^). 

s>0 s>0 

Fix g G iS(R'^) \ {0}. The adjoint operator of Vg, defined by (ld*-F, h) = {F, Vgh), 
can be written as 


( 20 ) 


V*F=f F{x,^)7i{x,^)gdxdf, 


The adjoint Vf maps the Banach space into M^'^(R'^), in particular it 

maps iS(R^^) into iS(R^) and the same for their dual spaces. If F = Vgf we obtain 
the inversion formula for the STFT 


( 21 ) 




(the same holds when replacing by iS(R“) or iS'(R“)). 


2.3. Phase functions and canonical transformations. Let a be as in the In¬ 
troduction, real-valued and satisfying ([2]). The related classical evolution, given by 
the linear Hamilton-Jacobi system, following our normalization can be written as 

'27rdtx{t,y,p) = -Vi:a{x{t,y,p),^{t,y,p)) 

27rdt^{t,y,g) = V^a{x{t,y,p),^{t,y,p)) 

x{0,y,v) = p, 

.^(0,2/,p) = p- 


(22) 
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The solution {x{t,y,ri),^{t,y,ri)) exists for every t G M. Indeed, setting u : = 
(x, ^), F{u) := {—'V^a{u),'Vxa{u)), the initial value problem fl2^ can be rephrased 
as 


(23) 


u'(t) ^ F(u{t)), u(fo) ^ Mo, 


in the particular case to = 0. Observe that our assumptions on a imply the bound¬ 
edness of d^^Fj, for every |q;|> 0, j = 2d, hence in particular F : —)■ R^'^ 

is a Lipschitz continuous mapping. The previous ODE is an autonomous ODE with 
a mapping F G C°°(R^'^ —?• R^'^) having at most linear growth, hence ||F(m)|| < 
1 -|- ||m||. Hence for each uq G R^'^ and to ^ ® there exists a unique classical global 
solution M : R —)■ R^'^ (in this case u G C°°(R —)■ R^'^) since F G C°°(R^'^ —)■ R^*^)) to 
The solution maps Stf,(t) : R^'^ C°°(R R^'^), defined by StQ(t)uo = u(t), 

and Stgito) = Id, the identity operator on R^'^, are Lipschitz continuous mappings, 
obeying St^it) = S'o(t — to) and the group laws 

(24) Soit)So{t') = So{t + t'), 5o(0) = Id. 

The mapping S'o(t) is a bi-Lipschitz diffeomorphism with SQ^{t) = S'o(—t). Follow¬ 
ing the notations of [TOl 120] , we call the bi-Lipschitz diffeomorphism 

(25) Xt{y,v) ■= So{t){y,v), {y,v) G R^'^. 

The theory of Hamilton-Jacobi allows to find a T > 0 such that for t g] — T, T[ 
there exists a phase function <h(t, x,//), solution of the eiconal equation 


(26) 


27idt^ + a(x, Va;<h) = 0 
<h(0, X, rj) = XT] 


The phase <h(t, x, rj) is real-valued since the symbol a(x, ^) is real-valued, moreover 
$ fulhlls the condition of non-degeneracy: 

(27) I det( 9 ^^^ 4 '(t,x,r 7 )| > c> 0, (t, x, 77 ) g] - T, T[xR^'^, 

after possibly shrinking T > 0, and satisfies 

(28) |afa"^<h(t,x,?7)| < Cfc,a, |q:| > 2 , fc > 0 , (t,x,77) g] - T,T[xR^'^. 

The relation between the phase $ and the canonical transformation % is given by 


(29) 


{x, VrHt, x,ri)) = Xi(V,#(f, X, t €] - T, T\. 


In particular. 


y{t,x,7]) = 
i{t,x,ri) = 


and there exists J > 0 such that 

dx 


Vr,<h(t,x,r7) 

Vx^{t,x,ri), 
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Observe that each component of Xt is a function with bounded smooth derivatives 
of any order in ] — T, Using we observe that the same holds in fact 

for every t G M. 

For f g] — T, T[, the phase function •) is a tame phase, and similarly for the 
canonical transformation Xt^ according to the following dehnition [TOl Dehnition 
2 . 1 ]: 

A real and smooth phase function on is called tame if: 

(i) For z = (x, 7]), 

(32) \df^{z)\<C^, |a|>2; 

(ii) There exists c > 0 such that 

(33) |det aj^<h(a;,77)| > c. 

The mapping defined by {x,f) = xiVihii which solves the system 

f y(x,?]) = Vr,^(x,?]) 

(34) 

[ = V^^(x,?]), 

is called tame canonical transformation. 

Observe that we have no assumption of homogeneity for large {x, rj), nevertheless 
the mapping x is well-dehned by the global inverse function theorem. The mapping 
X is a smooth bi-Lipschitz canonical transformation (i.e. it preserves the symplectic 
form) and satishes, for {x,f) = xidJ^v)) 


(36) 


ld;xi(z)\ + jal >\, z = 


The mapping x enjoys 


(36) \det —{y,r])\>6 

(that is (l3T|) for the canonical transformations of the Hamilton-Jacobi theory), 
which allows to uniquely determine, up to a constant, the related tame phase 
function (see [101 Section 2]). 


3. Fourier Integral Operators 

3.1. The classes FIO{x)- The class FIO{x) was introduced in (TU] and its deh¬ 
nition can be rephrased as follows. 

Definition 3.1. Let g G iS(M'^) be a non-zero window function. Consider a canon¬ 
ical transformation x which is a smooth bi-Lipschitz diffeomorphism and satisfies 
(l35|) . A continuous linear operator T : iS(M'^) —?■ S' (M'^) is in the class FIO{x) if 
its (continuous) Gabor matrix satisfies for all s > 0 the decay condition 

(37) \{Tn{w)g,r:{z)g)\<C{z-xiw))~f V2:,wGM^'^, 
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for a constant C > 0 depending on s. 

Note that we do not require fl36|) to be valid. 

The class F 10(E) = [J^ FIO{x) is the union of these classes where x runs over 
the set of all smooth bi-Lipschitz canonical transformations satisfying 0351) . 

The following properties are proved in m- 

(i) Boundedness ofT on ([101 Theorem 3.4]): 

If T G FIO{x), then T can be extended to a bounded operator on 
(in particular on L^(R'^)). 

(ii) The algebra property f [TI]l Theorem 3.6]): For i = 1,2, 

(38) T« G FIO{xi) ^ g FIO{xi ° X 2 ) • 

(hi) The Wiener property (Uni Theorem 3.7]): IfT G FIO(x) andT is invertible 
on L2(M^), then G FIO{x-^). 

These properties imply that the union FIOiE) is a Wiener subalgebra of £(L^(R‘^)), 
the class of linear bounded operators on L^(R‘^). Property (ii) can be rehned as 
follows. 

Lemma 3.2. For g FIO{xi), i = 1,2, the continuous Gabor matrix of the 
composition is controlled for every s >0 by 

(39) \{T^^''T^‘^\{w)g,Ti{z)g)\ < C'oCiC2(^ - Xi o X2(w))"^ w, 2 ; G R^'^, 

where Ci > 0 is the constant of in dSZD, ^ = 1,2, whereas (Pq > 0 depends only 
on s and on the Lipschitz constants of xi Xi^ ■ 

The proof is an easy consequence of Lemma 2.11 in [20] . 

By induction we immediately obtain 

Proposition 3.3. For n E N, n > 2, g FIO{xi)j i = 1, ■ ■ ■ , n, we have 

(40) • • •T("V(w)5(,7r(^)5()| < CqCi ■ --Cn^z - yi o ya o • • • o Xn{w))~\ 
where Cq depends on s and on the Lipschitz constants of the mappings: 

Xl, X^^ Xl ° X2, (Xl ° X2)~\ • • • , Xl O X 2 O • • • O Xn-U (Xl O X 2 O • • • O Xn-l)”^ 

Observe that, using Schur’s test and the same techniques as in the proof m 
Theorem 3.4], it is straightforward to obtain the following weighted version of [TU] 
Theorem 3.4]. 

Theorem 3.4. Consider pi G Then for every 1 < p < 00 , T G FIO{x) 

extends to a continuous operator from into MJ^. 
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Let us underline that jj. o x E since Vr o x ^ '^r, due to the bi-Lipschitz 

property of x- In particular 

(41) T : ^ 1 < p < cx), r G R. 

If X = Id, then the corresponding Fourier integral operators are simply pseudodif¬ 
ferential operators, as already shown in |27j . 

The characterization below, written for pseudo differential operators in the Kohn- 
Nirenberg form a{x,D), works for any r-form (in particular Weyl form a'^{x,D)) 
of a pseudodifferential operator. 

Proposition 3.5. Fix g G iS(R'^) and let a G iS'(R^'^). For s G R, the symbol 
a belongs to only if for every s > 0, and for suitable constants 

C = Cs depending on s, 

(42) \{a{x,D)'K{w)g,'K{z)g)\<Cs{z-w)~\ 

Similarly, under additional assumptions on the classes FIO{x), their operators 
can be written in the following integral form {FIOs of type I): 

(43) J(a,«F)/(a;)= / p)/(p) dp, / G 5(R‘'), 

jR'i 

where a G S'gQ(R^'^)) and $ a tame phase function. More precisely, this particular 
form is allowed starting from the class FIO{x) whenever the mapping x enjoys the 
additional property (l36)) as explained in the following characterization [TOl Theorem 

4.3] . 

Theorem 3.6. Considerg G iS(R'^). Let I be a continuous linear operator S(Mf) —)■ 
iS'(R‘^) and x be a tame canonical transformation satisfying fl36|) . Then the follow¬ 
ing properties are eguivalent. 

(i) I = /(cr, <F^) is a FIO of type I for some a G 

(ii) /gf/0(x). 

For X = Id we recapture the characterization for pseudodifferential operators of 
Proposition 13.51 

Remark 3.7. We shall apply the preceding results to the mappings Xt{x, g) coming 
from the Hamilton-Jacobi system fl22l) . so we need to be more precise on the estimate 
dslD, namely we have to show how the constants C depend on the time variable t. 
To this end, in short: gluing together the results |9l Theorem 3.3] and [101 Theorem 

4.3] , we obtain that the constants C = Cg in fl37|l . s > 0, can be estimated in terms 
of the Sq Q-semi-norms of a. Provided the Lipschitz constants of Xt are uniformly 
bounded with respect to t, as we have in fl2^ . we may assume Cs{t) G C(] — T,T[) 
in (imi if the semi-norms of at can be estimated by functions in C(] — T,T[). 
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4. Proof of Theorems 11.11 and 11.21 

Let us consider the Cauchy problem ([T]) with a(z) real-valued satisfying (j2]). 
From the results recalled in the Introduction we have for the representation 
(jH). More precisely: 

Proposition 4.1. There exists a constant T > 0, a symbol a{t,x,ri) G C°°(] — 
T, r[, S'qq) a real-valued phase function $ satisfying and (E7|) such that the 
evolution operator can be written as 

(44) {e"^^uo){t,x) = {FtUo){t,x) 
where Ft is the FIO of type I 

(45) {FtUo){t,x)= [ e^"'*’(*’^’^V(t,x,r7)tro(h)c^^- 


Remark 4.2. Notice that the function *h(t, •) of Proposition f.l and the related 
canonical transformation Xt in are tame. The Lipschitz constants of Xt nnd 


-1 


Xt " can be controlled by a continuous function oft on the interval ] — T,T[ and 
thus can be chosen uniform with respect to t on] — T,T[. 

From Theorem 13.61 we infer that the propagator belongs to FIO{F) for every 

hxed t e] — T,T[. More precisely: 


Proposition 4.3. 
(46) 


Under the assumptions of Proposition \4^ we have 

e"«eF/0(x,). tel-r,r[. 


where Xt is defined in (12B]) . Moreover for every s > 0 there exists C(t) 
C(] — T,T[) such that, for every g G iS(M'^) the Gabor matrix satisfies 

(47) \{e^^^7r{w)g,7r{z)g)\<C{t){z-Xt{w))~f 


Cs{t) G 


The last part of the statement follows from Remark 14.21 and Remark 13.71 
The previous proposition gives a representation of for |t| < T. Using the 
group property of the propagator we may obtain an expression of e®*^ for every 
t G M. Indeed, a classical trick, jointly with the group property of e®*'^, applies. 
Namely, we consider Tq < T/2 and dehne 


Ih — ]hTo, [h -|- 2)To[, h G Z 
For f G J/i, by the group property of e®*^: 

(48) 


^itH _ ^i{t-hTo 


)H (^^i{hTo)H/\h\yh\ 


and using Proposition 14.11 one can write 


itH 


\h\ 
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In general, or even the composition cannot be represented as 

|h| 

a type I FIO in the form (jl]). We shall prove below that the evolution is in the 
class FIO{xt) for every t G M, with y dehned in fl 2 ^ . so that this class is proven 
to be the right framework for describing the evolution 

Theorem 4.4. Given the Cauchy problem o with a{z) real valued satisfying ([2|), 
consider the mapping Xt defined in fl25|) . Then 

(50) e FlOixt), t G M 
and for every s > 0 there exists C(t) G C(M) such that 

(51) \{e^^^Tr{w)g,TT{z)g)\<C{t){z-Xtiw))~fi w, z G t G M. 

Proof. We fix To < T/2 as above. For f G M, there exists h & "L such that t ^ Ih- 
Using Proposition 14.31 for p = t — hTo g] — T,T[ we have that G FIO{xti) 
and for t 2 = jj^To g] — T,T[, G FIO{xt2)^ and for every s > 0, there exists 
a continuous function C{t) on ] — T,T[ such that (1471) is satisfied for t = ti and 
t = t 2 . Using the algebra property fl38l) . we have 

^ pio{xt, o (xij'"') 

and the group law (|M)l for Xt{y,v) = So(t){y,ri) gives 

Ail ° {xt2)^ ^ Xh + \h\t2 Xty 

as expected. Then, using fl40|) we obtain that the Gabor matrix of the product 
ptiH^^it 2 Hyh\ jg controlled by a continuous function Ch{t) on Ih- Finally, from the 
estimates 

\{e^^^7r{w)g,7r{z)g)\ < Ch{t){z - Xt{w))~fi t G 4, 

with Ch G C{Ih), it is easy to construct a new continuous controlling function C(t) 
on M such that (|5T|l is satisfied. □ 

In particular, the estimate fl5Tl) gives Theorem ll.il Using Theorem 13.41 we obtain 
Theorem 11.21 

5. Gabor singularities and proof of Theorem 11.31 

We want now to localize in the Gabor singularities of a distribution and 
study the action on them of 

For F C we define the h-neighborhood F^, 0 < 5 < 1, as 

(52) F 5 = { 2 ; G R^'^ : \z — Zq\ < 6{zo) for some Zq G F}. 

We begin to list some properties of the ^-neighborhoods, for the proofs we refer to 
[I 8 l Lemmas 7.1, 7.2]. 
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Lemma 5.1. Given 6, we can find 6*, 0 < 5* < 5, such that for every T C 
(53) cr^, 


(54) 

Lemma 5.2. Let x be a smooth bi-Lipschitz canonical transformation as in the 
preceding sections. For every 6 there exists 5*, 0 < 5* < 6, such that for every 

r c 

(55) x(r5*) c x(r)5, 

(56) x(r)<5* c x(r5). 

The constant S* depends on x and 6 but it is independent ofT. 

In the following we shall argue on / G iS'(M'^), and take windows g G iS(R'^). 
Since Lls>oM^^{R‘^) = we have for some Sq > 0 

(57) [ \VJ{z)\P{z)-P^°dz<oo. 

Definition 5.3. Let f G iS'(M'^), g G iS(M'^) \ {0}, T C 1 < p < oo, r G M. 
We say that f is Mf -regular in V if there exists 5 > 0 such that 

(58) [ \Vgf{z)\P{z)P'’dz < oo. 

JTs 

(obvious changes if p = oo). 

Of course, fl55D gives us some nontrivial information about / only when F is 
unbounded. We shall prove later that Definition 15.31 does not depend on the choice 
of the window g G iS(M'^). 

Theorem 5.4. Let and Xt be defined as in the previous Sections, fix Uq G 
iS'(M‘^) and T C If Uq is Mf-regular in F, then e^^^Uo is Mf-regular in Xt(F). 

Proof. For sufficiently small 5 > 0 we have fl55D in F^ whereas dSZD is valid in 
for some Sq > 0. Now, from Theorem II.II we have 

(59) Vg{e^^^Uo){z) = j k{t,w, z)VgUo{w) dw 
with 

(60) \k{t,w,z)\ = \ {e^^^n{w)g,7r{z)g)\ < {z - Xt{w))~fi 
for every s > 0. 
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We want to show that is regular in Xi(r). To this end, using (]56|1 in 

Lemma [521 we take hrst 6* < 6 such that Xt(r)5* C Xt(r5) and then using fl55]) in 
Lemma [5.1! we £x S' < S* such that 

(61) (xt(r)y)^, c Xi(r)5* c xi(r5). 

Note that for w ^ L^, i.e. Xt{w) ^ Xi(r5), and z G Xi(r)5' we have 

(62) Iz - Xi(w)| > max{(z), (w)} 

since Xt{u!) ^ ((Xi(r)5')^, in view of flMD . and we may use as well fISTD . 

Assuming for simplicity p < oo, we shall prove 

(63) / \Vg{e'^^uo)iz)\P{zrdz<oo, 
with S' determined as before. Using (15^ and (160]) . we estimate 

(64) \{zyVg{e'^^Uo){z)\< [ I{z,w)dw, 

JR2d 

with 


(65) I{z,w) = {zY{z - xt{w)) ^\VgUo{w)\. 

To show that is M^-regular in Xt(r) it will be sufficient to show that 


/(•, w) dw 


< oo. 




First, we estimate f^ 2 d I{z, w) dw for z G At (0)5' • We split the domain of integration 
into two domains 05 and \ 05. In \ 05 we use fl62l) to obtain 

[ I{z,w)dw<[ {zY{wy°{z - dw 

jR^-ixrs {wy° 


< 

n>j 


f {z -x,{w)r'‘- yyy dw 

R^'i\Ts {w) ^ 

\vM-)\' 




(•) 


.\so 


(^). 


So by (1071) and choosing s in (1001) so that r + sq — s < —2d, 


/R2d\r5 


/(•, w) dw 


< 11^.^'■+*0 ®||^i(-g 2 d)|| |V^Mo|(') '^“||LP(R2d) < OO- 


LPiXtiDs,) 















GABOR ANALYSIS FOR SCHRODINGER EQUATIONS 


15 


In the domain we have 


f I{z,w)dw< [ {zY{w) "\VgUo{w)\{wY dw 

Ts JTs 

^ / {z - Xt{w)Y~"\VgUo{w)\{wY dw 
drs 

<(xr'(-)r**(Charr,-|V^,«o|(-n (^) 


where Charp^ is the characteristic function of the set F^. The assumption fl58D 
yields to the estimate 



/(•, w) dw 


~ WiXt ^{■)Y 1lLi(R2d)|| \VgUo\{-Y\\LP{rs) 

Lnxt(r)s^) 


- IKT 1lLi(R2'i)||^3Mo(-)'’IUp(r5) < oo, 


because Xt is a bi-Lipschitz diffeomorphism and we may take s in flbOp so large that 
r — s < —2d. This concludes the proof. □ 


Remark 5.5. It is now easy to check that Definition \5.3\ does not depend on the 
choice of the window g G In fact, assume that the estimate (j58jl in Defini- 

tion \5.!^ is satisfied for some <5 > 0, and some choice of g E Then (j58|l is 

still satisfied, for some new 6 > 0, if we replace g with h G To prove this 

claim, observe that for every s > 0, 

\VHg{z)\ < {z)-Y zeR^Y 

Lemma l2.1\ then gives 

|i4/(w)| < J{w- z)-YVgf{z)\dz. 


The claim follows by splitting the domain of integration into F^ and \ F^, and 
arguing as in the proof of the preceding Theorem\5.4\ being now Xt = identity. 


The following definition allows one to describe the position in phase space of the 
singularities of a function /. 


Definition 5.6. Given f E S'(RY, we shall call filter of the Gabor singularities of 
f the collection of subsets 

^P(/) = {AcR‘^^: f IS Mf-regular mV = \ A}, 

cf. Deftnition \5.tA 


Let us write d{f) = drif) short. 

d{f) is a filter since if A G d{f) and A C A', then also A' G d{f), and moreover if 
Ai,..., A„ G d{f) then also E d{f)- Note that any neighborhood of oo i.e. 
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the complementary of a bounded set, belongs to d{f)- We have / G if and 

only if 0 G that is equivalent to saying that there exists Ai, ... ,An G 

such that = 0. 

Theorem 15.41 now follows. Indeed, the inclusion Xt{d{f)) C is just a 

restatement of Theorem 15.41 The opposite inclusion is equivalent to f) ^ 

d{f), namely to xT^did) C which is true by reversing the time. 
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